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Abstract 

Starting from the canonical formalism of relativistic (timeless) quantum mechanics, the formula- 
tion of timeless path integral is rigorously derived. The transition amplitude is reformulated as the 
sum, or functional integral, over all possible paths in the constraint surface specified by the (rela- 
tivistic) Hamiltonian constraint, and each path contributes with a phase identical to the classical 
action divided by K. The timeless path integral manifests the timeless feature as it is completely 
independent of the parametrization for paths. For the special case that the Hamiltonian constraint 
is a quadratic polynomial in momenta, the transition amplitude admits the timeless Feynman's 
path integral over the (relativistic) configuration space. 
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I. INTRODUCTION 



The idea that quantum mechanics can be well defined even if the notion of time is absent 

nn nn 

has been proposed |2] and developed in a number of different strategies [3H7J . The motiva- 
tion for formulating quantum mechanics in timeless description comes from the research on 
quantum gravity, as in the quantum theory of general relativity, the spacetime background 
is not fixed and generally it is not possible to make sense of quantum variables "at a moment 
of time" . This is closely related to the "problem of time" in quantum gravity [s| . 



ess) quantum mechan- 
9]. The formulation is 



In particular, a comprehensive formulation for the relativistic (time 
ics and its probabilistic interpretation are presented in Chapter 5 of 
based on the canonical (Hilbert spaces and self-adjoint operators) formalism and we wonder 
whether it also admits the covariant (sum-over-histories) formalism. In the conventional 
nonrelativistic (with time) quantum mechanics, the transition amplitudes are the matrix 
elements of the unitary evolution generated by the Hamiltonian and can be reformulated as 
the sum over histories, called the path integral (see |l0| for a detailed derivation). In the 
relativistic quantum mechanics, however, the concept of time evolution is not well defined in 
the fundamental level; therefore, conceptual issues and technical subtleties arise when one 
tries to derive the timeless path integral from the canonical formalism. 

The aim of this paper is to rigorously derive the timeless path integral for relativistic quan- 
tum mechanics, starting from the canonical formulation in 9j. It turns out the transition 
amplitude can be reformulated as the sum, or functional integral, over all possible paths in 
the constraint surface S specified by the (relativistic) Hamiltonian constraint H(q a ,p a ) = 
for the configuration variables q a and their conjugate momenta p a , and each path contributes 
with a phase identical to the classical action divided by h. Unlike the conventional path 
integral in which every path is parameterized by the time variable t, the timeless path in- 
tegral is completely independent of the parametrization for paths, manifesting the timeless 
feature. Furthermore, for the special case that the Hamiltonian constraint is a quadratic 
polynomial in p a , the timeless path integral over £ reduces to the timeless Feynman's path 
integral over the (relativistic) configuration space. 

The timeless path integral for relativistic quantum mechanics is appealing both conceptu- 
ally and technically. Conceptually, timeless path integral offers an alternative interpretation 
of relativistic quantum fluctuations and is more intuitive than the canonical formalism for 
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many aspects. It can give a new point of view about how the conventional quantum me- 
chanics with time emerges within a certain approximation and thus may help to resolve the 
problem of time. Technically, timeless path integral provides tractable tools to compute (at 
least numerically or approximately) the transition amplitudes which otherwise remain formal 
in the canonical formalism. For example, the semiclassical approximation for the timeless 
path integral can be developed a la the Wentzel-Kramers-Brillouin (WKB) method. 

In the research of loop quantum gravity (LQG), the sum-over-histories formulation is 
an active research area that goes under the name "spin foam models" (SFMs) (see and 



references therein for LQG and SFMs). In particular, over t 
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re past years, SFMs in relation 



14j . However, the Hamiltonian 



to the kinematics of LQG have been clearly established 
dynamics of LQG is far from fully understood, and although well motivated, SFMs have not 
been systematically derived from any well-established theories of canonical quantum gravity. 



Meanwhile, loop quantum cosmology (LQC) has recently been cast in a sum-over- 



listories 



formulation, providing strong support for the general paradigm underlying SFMS 15 



16j. In 



this paper, the timeless path integral is systematically derived from the canonical formalism 
of relativistic quantum mechanics, and we hope it will shed new light on the issues of the 
interplay between LQG/LQC and SFMs. 

This paper is organized as follows. We begin with a review on the classical theory 
of relativistic mechanics in Sec. [TT] and then a review on the quantum theory of relativistic 
mechanics in Sec. lIIIl The main topic is presented in Sec. lIVl where the timeless path integral 
is derived and investigated in detail. Finally, conclusions and outlooks are summarized and 
discussed in Sec. [V] 



II. CLASSICAL THEORY OF RELATIVISTIC MECHANICS 

The conventional formulation of classical mechanics treats the time t on a special footing 
and therefore is not broad enough for general-relativistic systems, which treat time on the 
equal footing as other variables. To include general-relativistic systems, we need a a more 
general formulation with a new conceptual scheme. A timeless formulation for relativistic 
classical mechanics is proposed for this purpose and described in detail in Chapter 3 of Q], 
excerpts from which are presented in this section with some new materials added to give a 
review and define notations. 
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A. Hamiltonian formalism 



Let C be the relativistic configuration space coordinated by q a for a = 1, 2, • • • , d with q a 
being the partial observables and d being the dimension of C. In nonrelativistic mechanics, 
one of the partial observables can be singled out and treated specially as the time t, i.e. 
q a = (£,<?*), but this separation is generally not possible for general- relativistic systems. An 
observation yields a complete set of q a , which is called an event. In nonrelativistic mechanics, 
an observation is a reading of the time t together with other readings q l . 

Consider the cotangent space Q = T*C coordinated by q a and their momenta p a . The 
space fl carries a natural one-form 6 = p a dq a . Once the kinematics (i.e. the space C of the 
partial observables q a ) is known, the dynamics is fully determined by giving a constraint 
surface £ in the space fl The constraint surface £ is specified by H = with a function 
H : Q — > R fc . Denote 7 an unparameterized curve in Q (observables and momenta) and 7 
its projection to C (observables only). The physical motion is determined by the function 
H via the following 

Variational principle. A curve 7 in C is a physical motion connecting the 
events q± and , if 7 extremizes the action 

S[l}= [p a dq a (2.1) 

in the class of the curves 7 which satisfy 

H(q a ,p a ) = 0, (2.2) 

(i.e. 7 G E) and whose projection 7 to C connect and q^- 

If k = 1, H is a scalar function and called the Hamiltonian constraint. If > 1, there 
is gauge invariance and if is called the relativistic Hamiltonian. The pair (C, if) describes 
a relativistic dynamical system. All (relativistic and nonrelativistic) Hamiltonian systems 
can be formulated in this timeless formalism. 

By parameterizing the curve 7 with a parameter r, the action (12. ip reads as 

% a ,Pa,iV,] = J dr L a { T ) d ^- - N t {r)H\q a ,p a )y (2.3) 

where the constraint (12.21) has been implemented with the Lagrange multipliers A^(r). Vary- 
ing this action with respect to Vj(r), q a {r) and p a {r) yields the constraint equation(s) (12. 2p 
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together with the Hamilton equations: 

= N :i{ r) dH Y' V "\ (2.4a) 
dr dp a 

dr n 1 dq a y 1 

For k > 1, a motion is determined by a A;- dimensional surfaces in C and different choices of 
the k arbitrary functions Nj(r) determine different curves and parametrizations on the single 
surface that defines a motion. For k — 1, a motion is a 1-dimensional curve in C and different 
choices of N(t) correspond to different parametrizations for the same curve. Different solu- 
tions of q a {j) and p a {r) for different choices of Nj(r) are gauge-equivalent representations 
of the same motion and different choices of Nj(t) have no physical significance. 
Along the solution curve, the change rate of H with respect to r is given by 
dH* dq a dH i dpadW dH-* dH i dH-* dH i 



dr dr dq a dr dp a dp a dq a dq a dp a 

= Ni (2.5) 

To be consistent, the physical motion should remain on the constraint surface S. That is, 
dH/dr has to vanish along the curve. Therefore, we must have the condition 

{IP, H j }\^ = 0, abbreviated as {IP, H j } w (2.6) 

for all % and j. A function F(q a ,p a ) defined in a neighborhood of £ is called weakly zero if 
F\ s = (abbreviated as F w 0) and called strongly zero if 



| s = and 



(dF 


dF\ 


\dq a ' 


Op a) 



= 0, abbreviated as F ~ 0. (2.7) 

s 

It can be proven that F ~ implies F ~ fiH t for some functions fi{q a ,p a )- Consequently, 
we have 

{H\W}^r\{q a ,p a )H k . (2.8) 



The condition (12.61) ensures all constraints H l to be first class. (See 17| for more about 
constrained systems and the concept of first class constraints.) 



B. Nonrelativistic mechanics as a special case 

The conventional nonrelativistic mechanics can also be formulated in the timeless frame- 
work as a special case. For the nonrelativistic systems, the relativistic configuration space 
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has the structure C = MxCo, where Co is the conventional nonrelativistic configuration space; 
i.e., q a = (t, q % ) as one of the partial observables is identified as the time t. Correspondingly, 
the momenta read as p a = (pt,Pi) with p t being the conjugate momentum of t and Pi being 
the conjugate momenta of q l . The Hamiltonian constraint is given by 

H(t, q l ,p t , Pi ) = Pt + H (q l , Pl ,t), (2.9) 

where Ho(q\pi, t) is the conventional nonrelativistic Hamiltonian function. Given the Hamil- 
tonian constraint in the form of (12.91) . the Hamilton equations (12.41) lead to 

* =W(T) , £ = (2.10.) 



which read as 



and 



dp t dH 

-r = — 2.11 

dt dt K J 

dc?_ = dHo dph = _dHo (2 u) 

dt dpi dt dtf' 1 ' ' 

if particularly we use t to parameterize the curve of solutions. Furthermore, the constraint 
( 12. 2 p dictates p t = —H Q . Thus, the momentum p t is the negative of energy and it is a 
constant of motion if H = Ho(q l ,pi) has no explicit dependence on t. The equations in 
( 12.121) are precisely the conventional Hamilton equations for nonrelativistic mechanics. 

The Hamilton equations in (I2.12p form a system of first-order ordinary differential equa- 
tions. Given the initial condition q l (t Q ) = g and Pi(to) = Pio at the time to, the existence 
and uniqueness theorem for ordinary differential equations states that there exists a solution 
of ( I2.12p given by q l = q l (t) and Pi = Pi(t) for t £ R, and furthermore the solution is unique. 1 
As a consequence, q l and evolve as functions of t, and a physical motion is an open curve 
in C = R x Co, along which the observable t is monotonic. 

A dynamical system in which a particular partial observable can be singled out as t 
such that the Hamiltonian is separated as in the form of (12. 9p is called deparametrizable. 
For deparametrizable systems, the change of t is in accord with the ordinary notion of time, 



1 In order to apply the existence and uniqueness theorem, we assume dH$/dq\ dH /dpi, d 2 H /dq l 
d 2 H /dpi 2 and d 2 H /dq l dpj all continuous. 
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which does not turn around but grows monotonically along the physical motion. Generically, 
however, relativistic systems might be non-deparametrizable — no preferred observable can 
serve as the time such that other variables are described as functions of time along the 
physical motion. The classical theory predicts the physical motion as an unparameterized 
curve, which gives correlations between physical variables, not the way physical variables 
evolve with respect to a preferred time variable. In the next subsection, we will introduce 
the timeless double pendulum as an example to illustrate the timeless feature. 



C. Example: Timeless double pendulum 

Let us now introduce a genuinely timeless system as a simple model to illustrate the 
mechanics without time. This model was first introduced in (3), Q| and used repeatedly as 
an example in 

Consider a mechanical system with two partial observables, a and b, whose dynamics is 
specified by the relativistic Hamiltonian 

H(a,b,p a ,p b ) = -1 (pl+pl + a 2 + b 2 - 2E) (2.13) 

with a given constant E. The relativistic configuration space is C = M? coordinated by a and 
b, and the cotangent space Q = T*C is coordinated by (a,b,p a ,pb). The constraint surface 
£ is specified by H = 0; it is a 3-dimensional sphere of radius \/2E in O. 
In the N(t) = 1 gauge, the Hamilton equations f)2.4p give 



da alb alp a ..,.„ 

T=Pa, T=Pbi ~r = - Q > ~~l = - & > ( 2 - 14 ) 

GET GET GST GST 

and the Hamiltonian constraint (12. 2p gives 

a 2 + b 2 + pl+pl = 2E. (2.15) 

The general solution is given by 

a(r) = A a sm(r), b(r) = A b sin(T + /3), (2.16) 

where A a = \/2E sin a and A b = \/2E cos a, and a and (3 are constants. 

Therefore, physical motions are closed curves (ellipses) in C = 1R 2 . (Choosing different 
gauges for iV yields the same curve with different parametrizations.) This system is non- 
deparametrizable and does not admit a conventional Hamiltonian formulation, because, as 



discussed in Sec. Ill physical motions in C = 1 x Co for a nonrelativistic system are 
monotonic in t and thus cannot be closed curves. 



D. Lagrangian formalism 



Consider the special case that the relativistic Hamiltonian is given in the form: 



(2.17) 



where a a , (3 a and 7 a are constant coefficients, and V(q a ) is the potential which depends only 
on q a . This form is quite generic and many examples of interest belong to this category 
such as the relativistic particle (free or subject to an external potential), the timeless double 
pendulum (harmonic or anharmonic) and the nonrelativistic system as described by (12. 9p 
with H = J2iPl/^ m i + Vitfit)- The Hamilton equations (I2.4p yields 



da a 

-j- = N{2a aPa + P a q a + la ) 

GIT 

AT (r M dV 



(2.18a) 
(2.18b) 



Equation fl2.18al) gives the relation between the momenta p a and the "velocities" q a : = 
dq a /dr, through which the inverse Legendre transform recasts the action (12. 3ft in terms of 
the Lagrangian function: 



S[q a } q a ,N;r]= J drL(q a ,q a ,N) 



dr[T — 



^ - M ~ la 



NV{q a ) 



(2.19) 



Variation with respect to N yields 

5S_ _ dL_ _ 
6N ~ dN ~ 
x ^ 1 



Aa a 



N 



-f3 a q a -la 



E 



q 



2a n N 



N 



- (3 a q a ~ la 



-V 



\ a a a / 



(2.20) 



In this subsection, the repeated index a is not summed unless is explicitly used. 
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which is precisely the Hamiltonian constraint ( 12. 2p . On the other hand, variation with 
respect to q a gives the equation of motion as a second-order differential equation: 

5S dL d dL 



5q a dq a dr dq a 







(M-)=py + P ala -2a a ^ (2.21) 



Ndr \NdrJ ™ ° dq a 

which is equivalent to (I2.18p . 



III. QUANTUM THEORY OF RELATIVISTIC MECHANICS 

The timeless formulation for relativistic classical mechanics is reviewed in Sec. [TTJ Based 
on the Hamiltonian framework of the classical theory, the quantum theory of relativistic me- 
chanics can be formulated in canonical formalism. Unlike the conventional quantum theory, 
relativistic quantum mechanics does not describe evolution in time, but correlations between 
observables. The timeless formulation for relativistic quantum mechanics is described in de- 



tail in Chapter 5 of |9|, excerpts from which are presented in Sec. IIII Al and Sec. IIII Bl to 
give a review. Issues on the physical Hilbert space are detailed in Sec. IIII CI and the physical 
interpretation of quantum measurements and collapse is discussed in Sec. IIII Dl 



A. General scheme 



Let C be the relativistic configuration space for the classical theory as described in the 
Sec. Ill Al The corresponding quantum theory can be formulated timelessly in the following 
scheme: 

Kinematical states. Let S C /C C S' be the Gelfand triple defined over C with the measure 
d d q a = dq 1 dq 2 ■ ■ ■ dq d . 3 The kinematical states of a system are represented by vectors 
l^) G /C, and K. is called the kinematical Hilbert space. 

Partial observables. A partial observable is represented by a self-adjoint operator in K,. 
The simultaneous eigenstates \s) of a complete set of commuting partial observables 

3 That is, S is the space of the smooth functions f(q a ) on C with fast decrease, JC = L 2 [C, d d q a ] is a Hilbert 
space, and S' is formed by the tempered distributions over C. 
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are called quantum events. In particular, q a and p a are partial observables acting re- 
spectively as multiplicative and differential operators on ip(q a ); i.e., q a, il){q a ) = q a tp(q a ) 
and f) a ip{q a ) = —ihdip{q a )/dq a . Their eigenstates \q a ) (defined as q a \q a ) = q a \q a )) and 
\p a ) (defined as p a \Pa) = Pa\Pa)) are both quantum events. 

Dynamics. Dynamics is defined by a self-adjoint operator H in JC, called relativistic Hamil- 
tonian. The operator from S to S' schematically defined as 

dre~ lTil (3.1) 

is called the "projector" . 4 The matrix elements 

W(s,s') := (s\P\s') (3.2) 

are called transition amplitudes, which encode entire physics of the dynamics. 

Physical states. A physical state is a solution of the quantum Hamiltonian constraint 
equation: 

H\f) = 0, (3.3) 

which is the quantum counterpart of (12. 2ft . Given an arbitrary kinematical state 
\ip a ) G S, we can associate an element (^J G S', defined by its (linear) action on 
arbitrary states \ipp) G S as 

(**J^> = Jdr (e^Val^) = (V>a|P|^), (3.4) 

such that (^ a | is a physical state, namely, a solution to (|3.3p . The solution space is 
endowed with the Hermitian inner product: 

(*^|*^):=(*^|^>, (3-5) 

which is called the physical inner product. The Cauchy completion of the solution 
space with respect to the physical inner product (-|-) is called the physical Hilbert 
space and denoted as %. 
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The integration range depends on the system: It is over a compact space if the spectrum of H is discrete 
and over a noncompact space if the spectrum is continuous. The operator P is a projector in the precise 
sense only if zero is a part of the discrete spectrum of H. 
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Measurements and collapse. If the measurement corresponding to a partial observable A 
is performed, the outcome takes the value of one of the eigenvalues of A if the spectrum 
of A is discrete, or in a small spectral region (with uncertainty) if the spectrum is 
continuous. Measuring a complete set of partial observables A\ simultaneously is 
called a complete measurement at an "instance", 5 the outcome of which gives rise to 
a kinematical state \ip a ) (which is a simultaneous eigenstate of A; L if the spectra of 
Ai are discrete). The physical state is said to be collapsed to \^^ a ) by the complete 
measurement. 

Prediction in terms of probability. If at one instance a complete measurement yields 
\ip a ), the probability that at another instance another complete measurement yields 
\1pf3) is given by 

2 



V0c 



(3.6) 



where 



WtypM := (ii>{,m a ) = ds ds' M^) W(s,s')Ms')- 



(3.7) 



In particular, if the quantum events s make up a discrete spectrum, the probability of 
the quantum event s given the quantum event s' is 

2 



V ss > = 



W(s,s') 



(3.8) 



y/W{s',s') 

If the spectrum is continuous, the probability of a quantum event in a small spectral 
region R given a quantum event in a small spectral region R' is 

W(R, R') 



Vrr> 



y/W(R,R) y/W(R',R f ) 



where 



W(R,R') = [ ds I ds'W(s,s'). 

J R Jr> 



(3.9) 



(3.10) 



5 In the timeless language, a complete measurement is said to be conducted at some "instance", not at 
some "instant". 
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It should be noted that, unlike the classical theory, the relativistic quantum mechanics 
described above is not equivalent to the conventional quantum theory, even if the system is 
deparametrizable. In conventional quantum mechanics, the time t is treated as a parameter 
and not quantized as an operator. Thus, the measurement of t is presumed to have zero 
uncertainty (At = 0). In relativistic quantum mechanics, t is on the same footing as other 
observables q % and the measurement of t will yield nonzero At. For a simple harmonic 



oscillator govern 



was shown in 



jy the relativistic Hamiltonian H = p t + Hq = p t + p a /2m + mu a /2, it 



6) that, if At <C mAa 2 /h, we can ignore the temporal resolution At and 
idealize the measurement of t as instantaneous, and the conventional nonrelativistic quantum 
theory is recovered as a good approximation of the relativistic quantum mechanics. 

In the following, we will first take the timeless double pendulum as an example to see 
how the above scheme is carried out, and then elaborate on some intricate issues. 



B. Example: Timeless double pendulum 

Take the timeless double pendulum introduced in Sec. Ill CI as an example. The kine- 
matical Hilbert space is K = L 2 (M. 2 , dadb) , and the quantum Hamiltonian equation reads 
as 

H1>(a, b)= l - - h 2 ^ + a 2 -b 2 - 2E^j ^(a, b) = 0. (3.11) 

Since H = H a + Hb — E, where H a (resp. Hb) is the nonrelativistic Hamiltonian for a simple 
harmonic oscillator in the variable a (resp. b), this equation can be easily solved by using 
the basis that diagonalizes H a and Hb- Let 

^„(a) = (a|n) = -L= H n (a) e~ a2/2h (3.12) 
Vn! 

be the normalized nth eigenfunction for the harmonic oscillator with eigenvalue E n = h(n + 
1/2), where H n (a) is the nth Hermite polynomial. Clearly, the function 

*Pn a ,n b ( a i b ) ■= ^n a (a)^n 6 (&) = (a, b\n a , n b ) (3.13) 

solves fl3~TTT) if 

h(n a + n b + l) = E, (3.14) 
which implies the quantum theory exists only if E = h(N + 1) with iV e Z + U {0}. 



12 



Consequently, for a given N, the general solution of (13. lip is given by 

N 

V(a,b) = J2 c nMa)^N-n(b), (3.15) 

n=0 

and thus the physical Hilbert space H is an (N + l)-dimensional proper subspace of /C 
spanned by an orthonormal basis {\n,N — n)} n= o,...,jv- 

The projector P : S — > H is a true projector as W is a proper subspace of K for the case 
that the spectrum of H is discrete. Obviously, P is given by 

N 

P = ^2\n,N-n)(n,N-n\, (3.16) 

n=0 

which can be obtained (up to an irrelevant overall factor) from (13. lft : 

LI 



I ^ dre- lrfl cci/'drV \n a , n b ) e -ir(n a +n b+ i-E) ^ ^ | 

= ^ ^n +n 6 +i,js|n tt , n 6 ) (n a , n b \= P. (3.17) 



Here, the integration range is so chosen because exp(-irH) is periodic in r with period 
2n/h itE = h(N + l). 

The transition amplitudes are given by 

N 

W(a,b,a',b') := (a, b\P\a', V) = ^ (a, 6|n, N — n)(n, N — n\a, V) 

n=0 

^ -(a 2 +fe 2 +a' 2 +b' 2 )/2/l 

= E ^ \1 H n (a)H N ^ n (b)H n (a')H N _ n (b'), (3.18) 

' n!(iv — n ! 

n=0 v y 

which is the probability density of measuring (a, 6), given (a',b') measured at another in- 
stance. Furthermore, the probability of the quantum event (n a , nb) given the quantum event 
«, n' b ) is 

W[lpna,n b ,1pn' a ,n' b ] ■= (n a , Tl b \P\n' a , Tl' b ) = ^n.+ni^^^^. (3.19) 

C. More on the physical Hilbert space 

The operator P : S — >■ 5' maps an arbitrary element of 5 to its dual space 5". If zero is 
in the continuous spectrum of H, P maps S to a larger space S' and thus is not really a 
projector. In this case, the physical state (^^ a | mapped from \ip a ) is a tempered distribution. 
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P becomes a true projector only if zero is a part of the discrete spectrum of H such as in 
the timeless double pendulum. 



The construction in (13. ip is a special case for the group averaging procedure [18|, [19|, the 
idea of which is to averaging over all states along the gauge flow (generated by the constraint 
operator) to yield the physical solution which satisfies the constraint equation. In the special 
case, let \E) be the eigenstate of H with eigenvalue E, then schematically we have 

H\^ a ) = J drEe- lr{l \i) a ) = J dr J dE He- irA \E)(E\^ a ) 

= J drj dE E e~ irE \E) (E\ip Q ) = J dE 5(E) E\E) (E\ip a ) = 0, (3.20) 

thus showing that P maps an arbitrary kinematical state \ip a ) to a physical state which 
satisfies the constraint equation (I3.3p . Furthermore, it can be easily shown that (^JV^) = 
(^vJV^) if 1 = i^tp' I; an d therefore the physical inner product in H3. 51) is well defined. 

If there are multiple constraints, we have to solve the multiple constraint equations si- 
multaneously: 

iT|<0)=O, for i = 1,... ,k. (3.21) 
In the simplest case that [H\ W\ = for all i,j, the projector can be easily constructed via 



P = J dn--- J dr k e~ lntl (3.22) 

as a direct extension of f)3.ip . In general, however, H l do not commute, as classically the 
Poisson brackets vanish only weakly [see (12.61) and (12.81) ]. 

In the case that H l do not commute but form a closed Lie algebra, i.e., 

[H\H j ] = f\H k (3.23) 

with f l3 k being constants, the exponentials of H % form a Lie group G and the physical state 
can be obtained by group averaging: 



G 



l**J= / dn(U)U\ij; a ), (3.24) 
where dfi is the Haar measure. It follows 



c/'l^j = / dii(u)u'u\ii> a )= / dniir'-'u") u"\if> a ) 

Jg Jg 

dii(U")U"\ip a ) = \% a ) (3.25) 
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for any U' G G. The fact that |^ a ) is invariant under any U' G G implies that it is 
annihilated by the generators of G, namely, H 1 \ \l/,/,„ ) = 0. Furthermore, the physical inner 
product in (13 . 5 p is again well defined. (See [19| for more details and subtleties.) The 
averaging in ( 13 .22 j) is indeed a special case of (I3.24p . 

Generically, j % \ are functions of q a and p a in (I2.8p . and, correspondingly, H l do not form 
a closed Lie algebra in the kinematical space /C. In this case, it is much more difficult to 
obtain the physical solutions and to construct the quantum theory which is free of quantum 
anomalies (see [20( for the issues of anomalies). 



D. Remarks on measurements and collapse 

Imagine that a quantum system is measured by Alice and Bob at two different instances, 
yielding two outcomes corresponding to \ip a ) and |^), respectively. From the perspective of 
Alice, the physical state is collapsed to \^^ a ) by her measurement and Bob's measurement 
affirms her prediction. Bob, on the other hand, regards the physical state to be collapsed to 
l^ipa) by his measurement and predicts what Alice can measure. The striking puzzle arises: 
Who, Alice or Bob, causes the physical state to collapse in the first place? 

In the timeless framework, it turns out to be an invalid question to ask who collapses 
the physical state first, since we cannot make any sense of time. The seemingly puzzle 
is analogous to the Einstein-Podolsky-Rosen (EPR) paradox, in which a pair of entangled 
particles are measured separately by Alice and Bob. In the context of special relativity, if the 
two measurements are conducted at two spacetime events which are spacelikely separated, 
the time-ordering of the two events can flip under a Lorentz boost and thus has no physical 
significance. Alice and Bob can both claim that the entangled state is collapsed by her/his 
measurement and thus have different knowledge about what the physical state should be, yet 
the predictions by Alice and Bob are consistent to each other. In our case, the measurement 
at an instance is analogous to the measurement on a single particle of the EPR pair; the 
kinematical state is analogous to the (local) state of a single particle; and the physical state 
is analogous to the (global) entangled state of the EPR pair. A complete knowledge (usually 
from measurement) about the local state will collapse the global state at once through the 
entanglement, which is analogous to the dynamics (or say, transition amplitudes) in our 



case. Consistency also holds in our case as (Nt^J^) = fif^Jipa) ■ (See [21] for more on the 



15 



5PR paradox in the relational interpretation of quantum mechanics and also Section 5.6 of 
for more on the philosophical issues.) 



As a side remark, exploiting further the close analogy between the EPR pair and the 
timeless formalism of relativistic quantum mechanics, one might be able to conceive an 
analog of the Bell's inequality, which would help to elaborate on the interpretational and 
conceptual issues of relativistic quantum mechanics at the level of thought experiments. 

IV. TIMELESS PATH INTEGRAL 

The canonical formalism for relativistic quantum mechanics is described in Sec. IHIi All 
information of the quantum dynamics is encoded by the transition amplitudes (13. 2p . In 
particular, by choosing \s) = \q a ) and \s') = \q' a ), all physics can be obtained from the 
following transition amplitudes 



From now on, we will use the notation ~ to denote the equality up to an overall constant 
factor which has no physical significance, as any overall constant is canceled out in the 
numerator and denominator in (13. 8\\ . 

As a special case of group averaging, the integration range of r is taken to be a compact 



correspondingly integrate r over (—00,00). The unwrapping only gives rise to an overall 
multiplicative factor (which is divergent if not properly regulated). Therefore, in any case, 
up to an irrelevant overall factor, transition amplitudes can be computed by 



where (q a \e lTH \q' a ) can be thought as the transition amplitude for a kinematical state \q' a ) to 
"evolve' to the state \q a ) by the "parameter time" r. Equation (14. 2p sums over (q a \e tTH \q' a ) 
for all possible values of r, suggesting that W(q a , q' a ) is intrinsically timeless as the parameter 
time r has no essential physical significance. 




(4.1) 





(4.2) 
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Rigorously, the integration should be regularized via 



W(q a , q' a ) ~ lim 



yl M drl ^ a \ e ~ iTH \ c ( a ) 




(4.3) 



as a cut-off M is introduced to regulate the integral and the irrelevant overall factor to be 
finite. As we will see, the variable r corresponds to the parametrization of curves in the 
path integral and integrating over all r indicates that the parametrization of curves has no 
physical significance. The above regularization scheme is physically well justified, as it cuts 
off the curves in the path integral which are "too wild" (noncompact curves), given the 
endpoints q' a and q a fixed. 

In the following, starting from Q4.2j) . we will first derive the timeless path integral for the 
case of a single Hamiltonian constraint and then investigate it in more detail. In the end, 
we will study the path integral with multiple relativistic constraints. 

A. General structure 

For a given r, let us introduce a parametrization sequence: r = 0, Ti, t 2 , • • •, rjv-i, tn = r 
with Ti G R, and define Ar n := r n — r n _i. The conditions on the endpoints (to = and 
r n = t) correspond to ^2^=1 ^ r n = r - The mesh of the parameter sequence is defined to be 
max n= i j ... ^{lA-Tnl}. The parameter sequence is said to be fine enough if its mesh is smaller 
than a given small number e. 6 

As r is fixed now, identifying q a = q^ and q' a = q%, and using J2n=i = r ' we can 

6 Let X be a topological space and s£ [0,1]. A continuous map 7 : s 1— > j(s) <E X is called a path with 
an initial point s(0) = x$ and an end point s(l) = xi. The image of 7 is called a curve, which can be 
reparameterized with respect to a new variable r as 7 : t4 7(t) by introducing an arbitrary continuous 
function r : s n- r(s) G M. The parametrization sequence To = 0, t%, T2, ■ ■ •, tjv-i, tjv = T can be viewed 
as a discrete approximation for the reparametrization function r(s) with r(s = 0) = and r(s = 1) = r if 
we identify t„ = t(ti/N). For the case that r(s) is injective, the parametrization sequence is ordered (i.e. 
= To < Ti < ■ • • < tjv-i < tjv = t and Ar„ > if r > 0) and called a partition of the interval [0, r], which 
is used to define the Riemann integral as the continuous limit: JJ" /(r')<ir' = lim mcs h-s.o J^-Tp 1 /(t„)At„. 
In the timeless formulation of relativistic mechanics, a dynamical solution is an unparameterized curve 
in and its parametrization has no physical significance. In order to exploit the timeless feature, we 
should keep the parametrizations generic and not restrict ourselves to injective ones. Correspondingly, 
the partition is generalized to a parametrization sequence and the Riemann integral is generalized to the 
Riemann-Stieltjes integral as f Q /(s)dr(s) = lim mcs h->o ^2n=a /( n /A r )Ar„, which is well defined even if 
t(s) is not injective. 
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rewrite (q a \e lTH \q ,a } as 



(q a \e- irH \q' a ) = {q a N \e~ iATNH e~ i/lTN - lH ■ ■ ■ e~ iAriH \q^) 



'N-l 



_,a I — iA.TNH\ rl a 



q a N-i}(i a N-i\ e 



-iArjv— l-ff I „a 



<&- 2 > • • • <g?|e 



a I — iAtiH I a 



9&>, (4-4) 



Uj d d q° n] ^ m 
where we have inserted AT — 1 times the completeness relation 

J d d q a \q a ) (q a \ ■= J dq 1 ■ ■ ■ dq d {q 1 , ■ ■ ■ , q d ) (q\ ■ ■ ■ , q d | . (4.5) 

For a given arbitrary small number e, by increasing N, we can always make the parameter 
sequence fine enough such that meshjrj} < \r\/N < e. 7 Consequently, we can approximate 
each (qn\ e ~ lAl ~ nH \ ln-i) to the first order in e as 

(q a n+1 \e~ lAr ^\q a n ) = (q a n+l \l - iAr n+1 H(q a ,p a )\q a n ) + 0(e 2 ). (4.6) 

For the generic case that the Hamiltonian operator H is a polynomial of q a and p a and 
is Weyl ordered, with the use of the completeness relation for the momenta 



d d p a 

it can be shown that 



\Pa)(Pa | ■-- 



dpi--- dp d 



- - • ,Pd)(Pi,--- ,Pd\ 



(q a \H(q q ,p a )\q' a ) 



d d Pa 

(2nh) d 6XP 



h 



Pa{q a - q" 



H 



q a + q la 



,Pa 



(4.7) 



(See Exercise 11.2 in [l0| for the proof.) Applying (14. 8 p to (14. 6p . we have 

d d p r , 



a I — jAt„+ 



lH K) 



(2nh) 
d d p. 



e %w« +1 -«)A [i _ iAr n+1 /J((^ +1 + <£)/2,P»«] + 0(6 2 



e ip na Aq%/h e -iAr n+1 H{q^,p na ) _|_ (9( e 2 ) ; 



(4.9) 



(27T^) rf 

where we define q a n := + g£)/2 and Ag^ := g£ +1 - g£. 

Making the parametrization sequence finer and finer (by decreasing e or equivalently by 
increasing N) and at the end going to the limit e — > or A" — > oo, we can cast (14.41) as 

d d p n 

n=0 



'iV-1 



(g a |e-^|g ,a ) = Jim ( j d d q a n J ( J] 



v n=l 



'iV-l 



11 7 (2vr^) £ 

n=0 ^ v y 



JV-1 



exp - ^ Pna^l 



N-l 



x exp -i &T n+1 H(q°, 



Pr, 



(4.10) 



n=0 



7 More rigorously, for a given e, the large number N should be chosen to satisfy mesh-jV^} < |t|/JV < M/N < 
e, where M is the cut-off regulator defined in (|4.3|) . so that the C(e 2 ) term in (|4.9[) can be dropped for 
any value of |r|. In the end, we have to integrate (|4.10[) over all possible values of r to obtain W(q a , q' a ), 
and the regularization is essential to keep the 0(e 2 ) terms under control for arbitrary values of r. 
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In the limit N — > oo, the points q n and p n can be viewed as the sampled points of a 
continuous curve in Q = T*C given by j(t') = {q a {T'),p a {f')), which is parameterized by r' 
and with the endpoints projected to C fixed by q a (r' = 0) = q' a and q a (r' = r) = q a . That 
is, q n and p n are the sampled points of 7 as = q a {r n ) and p na = p a (r n ). In the treatment 
of functional integral, it is customary to introduce the special notations for path integrals: 

N-l 

II / d d q a n / Vq\ (4.11a) 



n=l 



^-1 » id 



Meanwhile, in the continuous limit (N — > 00), the finite sums appearing in the exponents 
in (I4.10p also converge to the integrals: 

and 

AT-l 

- 1 &T n+l H(p na , q a n ) -+ -i H(q a (r'),p a (r')) dr'. (4.13) 

n=0 ^ 

Note that the continuous limit above is defined via the Riemann-Stieltjes integral as an 
extension of the Riemann integral (see Footnote [6]). With the new notations, (14.101) can be 
written in a concise form: 

(q a \e- ir6 \q' a ) = J Vq a j Vp a exp J Pa dq a j exp \-i j H(q a (r'),p a (r')) dr'^j . (4.14) 

It is remarkable to note that (up to the factor i/h) the continuous limit in (I4.12p is 
simply the line integral of the one-form 9 = p a dq a over the curve 7, identical to (I2.ip . and 
is independent of the parametrization r. On the other hand, the integral in (14.131) depends 
on the parametrization of r. Thus, to compute W(q a ,q' a )) in (I4.2p . the integration over r 
only hits the second exponential in (14.141) and the first exponential simply factors out. The 
integration of the second exponential over r yields 



[°° drexp f-i [H(q a (r'),p a (r'))dr' 
J —00 V J 7 



<1tvx]> I ^zr l_H(q a (f),p a (f))dfj = 5 (J H (q a (f) , p a (f)) df ) , ( I7L0) 
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where we have rescaled the parameterization r' to f = t'/t so that the endpoints now read as 
q' a = q a (f = 0) and q a = q a (f = l). 8 The appearance of the Dirac delta function indicates 
that only the paths which satisfy f~ H(f)df = will contribute to the path integral for 
W(q a ,q' a ). The condition f H(f)df = is, however, still not geometrical, since f can be 
further reparameterized to f' = f'(f) to yield f~ H(f')df' ^ even with the initial and final 
values fixed, i.e., f'(f = 0) = and f'(f = 1) = 1. On the other hand, W(q a ,q' a ) cast 
in (14. 2p has no dependence on the parametrization whatsoever, which implies that, in the 
continuous limit, the contribution of a path 7 satisfying the condition f~ H(f)df = for a 
specific (rescaled) parametrization f is somehow exactly canceled by that of another path 
satisfying the same condition. In the end, only the paths restricted to the constraint surface 
(i.e., 7 G S, or equivalently H(r') = for all r' along the path) contribute to the path 
integral for W(q a , q' a ). The constraint 7 G S is now geometrical. 

How the aforementioned cancelation takes place is obscure. To elucidate this point, we 
exploit the fact that W(q a , q' a ) is independent of the parametrization and play the trick by 
averaging over all possible parametrizations. That is, up to an overall factor of no physical 
significance, we can recast W(q a , q' a ) by summing over different parametrizations as follows: 

W(q a ,q' a )~ fdr J [VAr]^^ (<f|e iT V> 
~ J dr j [VAr]^ Tn=T J Vq a j Vp a exp Jp a dq a ) exp 
where the notation j [PAr]^ Ar =T is a shorthand for 

t/N pt/N pt/N p 

dA Tl / ' dAr 2 --- / ' dAr N -> / [PAr] EArn=T , (4.17) 

■t/N J-t/N J-t/N J 




which sums over all fine enough (namely, meshjrj} < \r\/N) parametrization sequences for 
a given r. It is easy to show that 

/oo pt/N pt/N pt/N poo 

dr / ' dA n / ' dAr 2 --- / ' dAr N = IT / dAr n+1 , (4.18) 
-00 J-t/N J-t/N J-t/N n=Q J -oo 



En=l &T n =T 



8 In the expression of (|4.15j) . we have removed the cut-off regulator (i.e. the limit M — > 00 has been taken). 
More rigorously, we have removed the regulator before the limit N —t 00 is taken. The Dirac delta function 
in (|4.15[> would have been a nascent delta function if the regulator had not been removed. 



20 



when the cut-off regufator M is removed (also see Footnoted]). Consequently, for a given 
arbitrary parametrization r', renaming the varying Ar n as Ar n = /T^A^At^, we can rewrite 
KWf as 



W(q a ,q' a ) (4.19) 
Vq a [ Vp a [ VN exp (~ [ p a dq a ) exp \- % - ^ Ar^N^HiCPnaj , , 



n=0 



where we introduce the notation 

N-l 



n 



dN n+1 -»■ 



PAT. 



(4.20) 



n=0 " _0 ° 

Again, in the continuous limit, the finite sum converges to the Riemann-Stieltjes integral: 



N-l 



n=0 



h 



N(r')H(q a (r f ),p a (r'))dT\ 



(4.21) 



and (I4.19P can be neatly written as the path integral: 



W(q a , q a ) ~ J Vq a J Vp a J VN exp 
= / Vq a / Vp a / VN exp 



(4.22) 



Integration over A^ can be carried out to obtain the delta functional: 



N-l 



PA>x P ( - J N(r')Hdr'j ~ = Km J] 8(H(<%,p an )), 



(4.23) 



n=0 



and thus the path integral (14.221) can be written in an alternative form as 



W(q a ,q" 



Vq a / Vp a 5[H] exp 



h 



p a dq a 



(4.24) 



where insertion of the delta functional S[H] confines the path to be in the constraint surface 
(i.e. 7 G S). Note that the phase in the exponent in (14.241) is identical to the classical action 
defined in (12. ip (divided by K) and that in (I4.22p is identical to the classical action in (12. 3p 
with k = 1. Therefore, each path in S contributes with a phase, which is the classical action 
divided by h~. 

The path integral formalism is intuitively appealing. It gives us an intuitive picture about 
the transition amplitudes: W(q a ,q' a ) is described as the sum, with the weight exp(iS/h) 
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(where S is the classical action of 7), over all arbitrary paths 7 which are restricted to S 
and whose projection 7 to C connect q' a and q a . None of q a is restricted to be monotonic 
along the paths, and in this sense the formulation is called timeless path integral. The 
parametrization for the paths has no physical significance as can be seen in the expression of 
f!4.24p . which is completely geometrical and independent of parametrizations. On the other 
hand, the continuum notation of (14.221) is really a schematic for the discretized version: 

N-1 „ N-1 „ , d N-1 „ 

w) ~ n n ^ n / ^ 

n=l J n=0 J v ' n=0 J 

x ex P (~l 2 A <+i N n+iH{q°,Pna) j (4.25a) 



n=l 17 n=0 17 v 7 n=0 ^ 

x ex P ( ~\ H N n+1 H(CPna) j , (4.25b) 



n=0 



where Ar^ in f)4.25ap is absorbed to N n in (14. 2 5b j) and this only results in an irrelevant 
overall factor. The expression of ( I4.25bl) is explicitly independent of parametrizations. 9 

The contributing paths in the path integral can be very "wild" — not necessarily smooth 
or even continuous. This calls into question whether the path integral can achieve conver- 
gence. We do not attempt to present a rigorous derivation here but refer to {22] for the 
legitimacy issues and subtleties of the path integral. 

Each path in E contributes with a different phase, and the contributions from the paths 
far away from the stationary solution essentially cancel one another through destructive 
interference. As a result, most contributions come from the paths close to the stationary 
solution. The stationary solution can be obtained by taking the functional variation on 
(I4.22p with respect to N, q a and p a , which yield the classical Hamiltonian constraint (12. 2p 
and the Hamilton equations (12. 4p . Therefore, the stationary solution is the classical solution 
and we thus have the approximation 

W(q a } q' a ) w e * Sm > ( 4 - 26 ) 



9 Perhaps, more appropriately, the "timeless path integral" should be renamed "timeless 'curve' integral" , 
as in the rigorous terminology, a curve is defined as the unparameterized image of a path, which is specified 
by a parameter. However, we keep the name of "path integral" to conform the conventional nomenclature. 
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where 7, are the classical solutions which connect q' a and q a and S is the action. 10 Based on 
the path integral formalism, a semiclassical theory can be developed in the vicinity of the 
classical solutions a la the WKB method. 



B. Deparametrizable systems as a special case 



If the Hamiltonian happens to be deparametrizable, the classical Hamiltonian is in the 
form of (12.91) . and the path integral (I4.24p reads as 



W(q a , q" 



Vt J Vq l J Vp t J V Pi 5[p t + H Q ] exp 
Vt / Vq l I Vpi exp 



h 



(p t dt + Pidq 1 ) 



dq % 

h.L \ p ^- Ho]dt 



(4.27a) 



n=l 



f N-l 



- in /««.)( n/^'Aiin 



,n=l 



v n=0 



(2nh) d - 1 



x exp 



JV-l 



n=0 



At r 



(4.27b) 



On the other hand, in the conventional nonrelativistic quantum mechanics, the transition 
amplitude is given by the conventional path integral (see Chapter 11 of [IQfl): 



G{q\t-q'\t') := (g^e-^-'V*) 
i 



Vq % \ Vpi exp 



h.l, \ Pi %- H ° ]dt 



(4.28a) 



'N-l 



,n=l 



- >- ( n / ) [ n 



Ui jJ rt 



,71=0 



(2nh) d - 1 



x exp 



h 



N-l 



n=0 



J2( Pm ^-H ( ( f n ,p m ))At 



At 



(4.28b) 



where the path 70 is in the cotangent space T*C and its projection 70 on Co has the endpoints 
fixed at q l (t') = q H and q l (t) = q\ 

Equations (14.271) looks almost the same as (I4.28P except for the functional integral J Vt. 
In (I4.27p . all paths given by q a (r) = (t(r),g*(r)) (parameterized by an arbitrary parameter 
r) are summed and the quantum fluctuation in t is also included; by contrast, in (I4.28p . all 



10 Generally, there could be multiple classical solutions connecting q' a and q a (as in the case of the timeless 
double pendulum), if the system is not deparametrizable. 
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paths are given by q l {t) as t is treated as a parameter and with no fluctuation at all. In other 
words, (14.281) sums over only the paths which are monotonic in t, whereas (14.271) sums over 
all possible paths whether t is monotonic or not. The difference is profound and shows that 
the relativistic quantum mechanics is not equivalent to the conventional quantum mechanics 
even if the system is deparametrizable, as already commented in the end of Sec. IHI Al 

However, for most systems, we have the good approximation (14.261) and only the paths 
in the vicinity of the classical solution are important. Meanwhile, as discussed in Sec. Ill Bl 
the classical solution for a deparametrizable system is always monotonic in t. Therefore, in 
(I4.27p . it is a good approximation to sum over only the paths which are not too deviated 
from the classical solution and thus monotonic in t. In this approximation, (14.271) reduces 
to the conventional path integral (14.281) as J T>t factors out as an irrelevant overall factor. 
Therefore, the conventional path integral, although not equivalent to, is a good approxi- 
mation for the timeless path integral. Further research is needed to investigate when the 
approximation remains good and when it fails; this is closely related to the question when 
the fluctuation (uncertainty) of t can be ignored and thus the measurement of it can be 
idealized as instantaneous (see ^l, H| for the case of a simple harmonic oscillator). 



C. Timeless Feynman's path integral 



Consider the special case that the classical Hamiltonian is given in the form of (I2.17P and 
the Hamiltonian operator is Weyl ordered. As the Hamiltonian is a quadratic polynomial in 
p a , the path integral over Vp a in (I4.22p can be integrated out. That is, in the expression: 11 



W(q a ,q" 



(4.29) 



iV-l 



Vq a / VN ] ] 



n=0 



d d Pn 

(2ixh) c 



exp 



Af-l 



$>naAg£ - /\r' n+1 N n+l H{q a n . 



Pr, 



n=0 \ a 



the integration over each p na can be explicitly carried out: 



dp na exp - [pna^ql ~ Ar^ +1 N n+1 (a a p 2 na + (3 a p na q^ + ^aPna)] 



oc 



cxp 



h 

i A< +1 iV ra+1 



-i 2' 



&< +1 N n+1 



(4.30) 



ii 



In this subsection, the repeated index a is not summed unless ^ a is explicitly used. 
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by the Gaussian integral f^dxe ax2+ ^ x = (tt / 'a) l/2 e^ ' /Aa . Noting that dN n+1/ 
2 dy/ Nn+i and introducing the shorthand notation: 



n / d^N^; 



WN 



(4.31) 



we then have 



W(q a , q" 



Vq a 



x exp 



WN 
1 ST ( Nn+1 

n=0 \ a 



(4.32) 



&q a n 



&< +1 N n+1 



~ Paq a n ~ la 



- N n+1 V(q a n ) 



which written in the continuous form reads as 
W(q a ,q' a )~ 



V^lvVNe^-ldA^ 



jj-P a q a -la 



- NV(q a ) , (4.33) 



where the "velocity" q a := dq a /dr' is the continuous limit of Aq®/ Ar n+1 . 

Therefore, in the special case that the Hamiltonian is a quadratic polynomial in p a , the 
transition amplitude admits a path integral formalism over the configuration space, whereby 
the functional integration over N is modified as J U\/N. This is called the configuration 
space path integral or Feynman's path integral. The configuration space path integral ( 14 . 3 3 j) 
sums over all arbitrary paths 7 G C whose endpoints are fixed at q' a and q a , and each path 
contributes with a phase, which is identical to the Lagrangian function as given in (I2.19p 
(divided by K). The functional variations on (I4.33P with respect to y/N and q a yield the 
classical Hamiltonian constraint and equation of motion as in (12.201) and (12. 211) . 12 This 
shows again that the stationary solution is the classical solution and thus ( I4.26P is a good 
approximation. 



D. Path integral with multiple constraints 

If there are multiple constraints and the constraint operators H l commute, the projector 
is given by (I3.22p and ( 14. 2 p can be directly generalized as 13 

/*oo HOC 

W(q a ,q' a )~ dr 1 --- dr k (g a |e-^=i ^V")- (4-34) 



12 Note that 5W/5\/N = 2\/N SW/SN. 

13 In this subsection, the repeated index i is not summed unless J2i i s explicitly used. 
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If each H l is a polynomial of q a and p a and Weyl ordered, the linear sum H' = ^\ t 1 H 1 is 
also a polynomial and Weyl ordered. Thus, by replacing r with 1 and H with H' in (14.141) . 
it can be shown 



q' a ) 



Vq a / Pp a exp 



h 



p a dq a 



xexp i-iY, fj % H\q a (f),p a (f))df\ , (4.35) 



where r is a parameter for the curve 7 with g /a = q a {r = 0) and q a = q a (r = 1). Redefining 
r l Af n as Ar^, we then have 



q' a ) 



Vq a / Vp a exp ( - / p a (ig 



x 



exp J_H\q a (T*),p a (T*))dA , (4.36) 



As in the case with a single constraint, the first exponential in (I4.36P is independent 
of parametrizations for the curve 7, and for the second exponential we can play the same 
trick by summing over different parametrizations to get rid of the seemingly dependence on 
parametrizations. Following the same steps in Sec. IIV A\ for each i, we have 

(4.37) 
(4.38) 



dr< I [PAr*] EATA=Tl (g a |e iTl "V a > 



Vq a / Vp a / VNi exp 



for a given arbitrary parametrization r'. After summed over [DAr l ]^ ATj=T , for each i, ( I4.34j) 
yields 



W(q a ,q'i) ~ / Vq a / Vp a \[ / VN t exp 



i=l 
k 



Vq a j Vp a \[5[W] exp % J L j_Padq a 



(4.39a) 
(4.39b) 



which is the direct generalization of (I4.22p and (14.24)) . In the path integral, each path in E 
contributes with a phase, which is the classical action given in (1 2 . 3 p divided by h. Functional 
variation on (|4.39a)) with respect to N iy q a and p a again yields the classical equations (12. 2j) 
and flZH) . 
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V. SUMMARY AND DISCUSSION 

Starting from the canonical formulation in the timeless path integral for relativistic 
quantum mechanics is rigorously derived. Given in ( I4.24p . the transition amplitude is for- 
mulated as the path integral over all possible paths in the constraint surface S (through the 
confinement by the delta functional ^[-ff]), and each path contributes with a phase identical 
to the classical action J~p a dq a divided by h. The alternative expression is given in (I4.22p . 
which is the functional integral over all possible paths in the cotangent space Q = T*C as 
well as over the Lagrange multiplier N. The timeless path integral manifests the timeless 
feature of relativistic quantum mechanics, as the parametrization for paths has no physical 
significance. For the special case that the Hamiltonian constraint H(q a ,p a ) is a quadratic 
polynomial in p a , the transition amplitude admits the timeless Feynman's path integral over 
the paths in the configuration space C, as given in (I4.33p . 

The formulation of timeless path integral is intuitively appealing and advantageous in 
many respects as it generalizes the action principle of relativistic classical mechanics by 
replacing the classical notion of a single trajectory with a sum over all possible paths. In 
particular, it is easy to argue that the classical solution contributes most to the transition 
amplitude and thus ( I4.26P is a good approximation for generic cases since the stationary 
solution is identical to the classical one. In the vicinity of the classical trajectory, the 
semiclassical approximation a la the WKB method can be developed. Furthermore, time- 
less path integral offers a new perspective to see how the conventional quantum mechanics 
emerges from relativistic quantum mechanics within a certain approximation (as discussed 
in Sec. HVBp and may provide new insight into the problem of time. 



The formulation of timeless path integral can be directly extended for the dynamical 
systems with multiple constraints as given in ( I4.39p . if the constraint operators H l commute. 
For the case that H l do not commute but form a closed Lie algebra, the projector is no loner 
given by H3 .221) but we have to invoke ( 13. 24ft to obtain the physical state, which leads to 



W(q a ,q' a )~ J dfi(9)(q a \e-^\q ,a ), (5.1) 

where Q % are coordinates of the Lie group G generated by H l . Starting from ( 15. ip and 
following the similar techniques used in this paper, one should be able to obtain the timeless 
path integral, but the measure of the functional integral Y\.%=\ I appearing in (|4.39aj) 
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would have to be nontrivially modified as 8 l play the same role of r l in (14.341) but now the 
nontrivial Haar measure dfi is involved and the nontrivial topology of G has to be taken 
into account. For the case that H l do not form a closed Lie algebra, it is not clear how 
to construct the quantum theory which is free of quantum anomalies even in the canonical 
formalism. The timeless path integral may instead provide a better conceptual framework 
to start with for the quantum theory. 

Throughout this paper, we have focused on simple mechanical systems, but not field the- 
ories. In Section 3.3 of [9|], the canonical treatment of classical field theories which maintains 
clear meaning in a general-relativistic context is presented as a direct generalization of the 
timeless formulation for relativistic classical mechanics (see also 23) and references therein), 



and the corresponding quantum field theory is formulated in Section 5.3 of js]. The timeless 
path integral for relativistic quantum mechanics derived in this paper should be extended 
for the quantum field theory described in [9J. We leave it for the future research. 

Furthermore, as the timeless path integral is systematically derived from the well- 
controlled canonical formulation of relativistic quantum mechanics, we expect it to provide 
new insight into the issues of the connection between LQG/LQC and SFMs. Extending 
timeless path integral to field theories will be particularly helpful. 
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